翻訳と辞書 Words near each other ・ "O" Is for Outlaw ・ "O"-Jung.Ban.Hap. ・ "Ode-to-Napoleon" hexachord ・ "Oh Yeah!" Live ・ "Our Contemporary" regional art exhibition (Leningrad, 1975) ・ "P" Is for Peril ・ "Pimpernel" Smith ・ "Polish death camp" controversy ・ "Pro knigi" ("About books") ・ "Prosopa" Greek Television Awards ・ "Pussy Cats" Starring the Walkmen ・ "Q" Is for Quarry ・ "R" Is for Ricochet ・ "R" The King (2016 film) ・ "Rags" Ragland ・ ! (album) ・ ! (disambiguation) ・ !! ・ !!! ・ !!! (album) ・ !!Destroy-Oh-Boy!! ・ !Action Pact! ・ !Arriba! La Pachanga ・ !Hero ・ !Hero (album) ・ !Kung language ・ !Oka Tokat ・ !PAUS3 ・ !T.O.O.H.! ・ !Women Art Revolution Dictionary Lists mini英和辞書 mini和英辞書 Webster 1913 Latin-English FOLDOC Wikipedia English ウィキペディア

翻訳と辞書　辞書検索 [ 開発暫定版 ] スポンサード リンク adaptive stepsize ： ウィキペディア英語版 adaptive stepsize Adaptive stepsize is a technique in numerical analysis used for many problems, but mainly for integration. It can be used for both normal integration (i.e. quadrature), or the process of solving an ordinary differential equation. This article focuses on the latter. For an explanation of adaptive stepsize in normal integration, see for example Romberg's method. As usual, an initial value problem is stated: : $y\text{'}(t)\; =\; f(t,y(t)),\; \backslash qquad\; y(a)=y\_a$ Here, it is made clear that ''y'' and ''f'' can be vectors, as they will be when dealing with a system of coupled differential equations. In the rest of the article, this fact will be implicit. Suppose we are interested in obtaining a solution at point $t=b$, given a function $f(t,y)$, an initial time point, $t=a$, and an initial solution $y\_a=y(a)$. Of course a numerical solution will generally have an error, so we assume $y\_b\; +\; \backslash epsilon\; =\; y(b)$, where $\backslash epsilon$ is the error. For simplicity, the following example uses the simplest integration method, the Euler method. Note that the Euler method is almost exclusively useful for educational purposes; in practice, higher-order (Runge-Kutta) methods are used due to their superior convergence and stability properties. Recall that the Euler method is derived from Taylor's theorem with the intermediate value theorem and the fact that $y\text{'}(t)=f(t,y)$: : $y(t+h)=y(t)+hf(t,y(t))+\backslash fracf\text{'}(\backslash eta,y(\backslash eta)),\; \backslash qquad\; t\; \backslash le\; \backslash eta\; \backslash le\; t+h$ Which leads to the Euler method: : $y\_^=y\_n+hf(t\_n,y\_n)$ And its local truncation error : $\backslash tau\_^=ch^2$ : $y\_^\; +\; \backslash tau\_^=y(t+h)$ We mark this solution and its error with a $(0)$. Since $c$ is not known to us in the general case (it depends on the derivatives of $f$), in order to say something useful about the error, a second solution should be created, using a stepsize that is smaller. For example half the original stepsize. Note that we have to apply Euler's method twice now, meaning we get two times the local error (in the worst case). Our new, and presumably more accurate solution is marked with a $(1)$. : $y\_\}=y\_n+\backslash fracf(t\_n,y\_n)$ : $y\_^=y\_\}+\backslash fracf(t\_\},y\_\})$ : $\backslash tau\_^=c\backslash left(\backslash frac\backslash right)^2+c\backslash left(\backslash frac\backslash right)^2=2c\backslash left(\backslash frac\backslash right)^2=\backslash fracch^2=\backslash frac\backslash tau\_^$ : $y\_^\; +\; \backslash tau\_^=y(t+h)$ Here, we assume error factor $c$ is constant over the interval $(t+h)$. In reality its rate of change is proportional to $y^(t)$. Subtracting solutions gives the error estimate: : $y\_^-y\_^\; =\; \backslash tau\_^$ This local error estimate is third order accurate. The local error estimate can be used to decide how stepsize $h$ should be modified to achieve the desired accuracy. For example, if a local tolerance of $tol$ is allowed, we could let h evolve like: : $h\; \backslash rightarrow\; 0.9\; \backslash times\; h\; \backslash times\; \backslash frac|\}$ The $0.9$ is a safety factor to ensure success on the next try. This should, in principle give an error of about $0.9\; \backslash times\; tol$ in the next try. If , we consider the step successful, and the error estimate is used to improve the solution: : $y\_^\; =\; y\_^\; +\; \backslash tau\_^$ This solution is actually third order accurate in the local scope (second order in the global scope), but since there is no error estimate for it, this doesn't help in reducing the number of steps. This technique is called Richardson extrapolation. Beginning with an initial stepsize of $h=b-a$, this theory facilitates our controllable integration of the ODE from point $a$ to $b$, using an optimal number of steps given a local error tolerance. Similar methods can be developed for higher order methods, such as the Runge-Kutta 4th order method. Also, a global error tolerance can be achieved by scaling the local error to global scope. However, you might end up with a stepsize that is prohibitively small, especially using this Euler based method. If you are interested in adaptive stepsize methods that use a so-called 'embedded' error estimate, see Fehlberg, Cash-Karp and Dormand-Prince. These methods are considered to be more computationally efficient, but have lower accuracy in their error estimates. == References == 抄文引用元・出典: フリー百科事典『 ウィキペディア（Wikipedia）』 ■ウィキペディアで「adaptive stepsize」の詳細全文を読む スポンサード リンク 翻訳と辞書 : 翻訳のためのインターネットリソース Copyright(C) kotoba.ne.jp 1997-2016. All Rights Reserved.